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Abstract: The generalized Tsallis statistics produces a 
distribution function appropriate to describe the interior solar 
plasma, thought as a stellar polytrope, showing a tail depleted 
respect to the Maxwell-Boltzmann distribution and reduces 
to zero at energies greater than about 20 k B T. The Tsallis 
statistics can theoretically support the distribution suggested 
in the past by Clayton and collaborators, which shows also a 
depleted tail, to explain the solar neutrino counting rate. 
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By imposing a variational principle one can conve- 
niently generalize both Boltzmann-Gibbs statistics and 
standard thermodynamics Jl],|2| . The lack of adequacy of 
the Boltzmann entropy is related to the breakdown of 
Boltzmann-Gibbs statistics for systems with long-range 
interactions ||,|[ such as, among others, the interior so- 
lar plasma, where the particles are exposed to many-body 
collisions and the energy available to a particular pair of 
particles is not defined ||. 

Even for a non perfect plasma, the ideal gas approxi- 
mation is usually considered a good one. The Maxwell- 
Boltzmann (MB) distribution is believed to be highly 
correct and the many-body physics involved does not 
determine sensible deviations from the standard statis- 
tics used ||. However, non-Maxewellian (flattened) dis- 
tributions in plasmas heated by inverse bremsstrahlung 
(collisional) absorption of sufficiently strong electromag- 
netic fields have been predicted and recently measured 
; in D - 3 He fusion plasma experiments, a nonther- 
mal component of the ion distribution and the possible 
consequences on the nuclear rates have been investigated 
@0. 

Among the distributions that deviate from the MB 
one, we pose our attention on distributions coming from 
generalized statistics. We give a very brief outline of the 
Tsallis generalization of thermodynamics and statistical 
physics , suitable for describing systems with long- 
range interactions |||| . 

If we have a system with W microscopic states, each with 
probability fi > normalized as 



£/* = ! 



then the entropy is given by 



w 

i=l 



(i) 



(2) 



*e-mail: Quarati@polito.it 



where k and q are constants. In the limit as q approaches 
unity, the well known expression Si — —k^^^filogfi is 
recovered and we can fix k — k B (the Boltzmann con- 
stant), i.e. the Tsallis statistics reduces to the standard 
one as q — > 1. 

We have shown a generalization of the Tsallis statistics 
by using a kinetic approach, based on the Fokker-Planck 
equation, and we have given to the Tsallis parameter q 
the meaning of a measure of the deviation from constant 
behavior of the diffusion coefficient D(v), a quadratic 
function of the velocity v [p^|-p^[. The Tsallis gener- 
alized statistics is actually widely used in many different 
physical problems (we send the reader to Ref.s Jl5 16 
where many references on the different applications are 
quoted) . 

One of the first problems where Tsallis statistics has 
been applied is that of stellar polytropes . Nobre and 
Tsallis Q and other authors have recently discussed the 
physical needs for departure from Boltzmann-Gibbs sta- 
tistical mechanics and thermodynamics for gravitational 
like systems JlTj and anomalous diffusion Jl6],[L8| (this 
last subject could be of great interest in the study of so- 
lar core, due to the well known problem of the diffusion 
of light elements like Li and Be ||l9|| ). 

Plastino and Plastino || have sought help from Tsallis 
entropy to find sensible distribution functions for stellar 
polytropes while that of MB gives unphysical distribu- 
tion functions. 

They found a range of variability of q from a relation be- 
tween the polytropic index n and q deduced comparing 
two different but equivalent expressions of the distribu- 
tion. We obtain a different relation between n and q 
because we impose a special constraint on the solar grav- 
itational potential. 

The internal structure of the sun can be considered 
polytropic with index n ranging between the value 3/2 
(7=5/3) and 5 (7=6/5) || (7 is the adiabatic parame- 
ter). In the interior regions where hydrogen ionization 
is changing rapidly, 7 is taken to be very close to unity 
(n = 00, q — 1), therefore in these circumstances the MB 
distribution holds. 
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We want to apply generalized statistics to derive a dis- 
tribution function for the interior solar plasma, relevant 
in calculating the nuclear fusion reaction rates responsi- 
ble for the neutrino flux emitted by the sun |2(],|2l[] . The 
following arguments support the program of this work. 

It is well known that one of the main problems with 
solar physics is related to the detected neutrinos arriving 
from the sun. All the different experiments have con- 
firmed a deficit in the flux relative to the predictions of 
standard theory of nuclear physics [jl9|,^l]]. In particular, 
the neutrino flux from the reactions involving 7 Be and 
8 B, mainly due to collisions at energies higher than the 
effective energy of pp reactions (4.58 k B T), is much lower 
than that predicted by the standard models, which uses 
MB distributions @. 

In the recent past we have tried to contribute to the 
solution of the solar neutrino problem deriving, as steady 
state solutions of the Fokker-Planck equation, statistical 
distributions that differ from the MB distribution for a 
depleted tail at high energies p^JT^| . The distribution 
/ is also the solution of a Boltzmann equation without 
collisions (at equilibrium the distributions obtained with 
and without the collision term coincide). 
Non-Maxwellian distribution can give thermonuclear re- 
action rates r =< av > smaller than the standard ones 
p| , allowing a reduced flux of neutrinos from the sun 
interior, particularly at energies above few k B T. 

Very recently we have shown that the Tsallis distribu- 
tion is a particular case of the family of distributions we 
derived by means of the Fokker-Planck equation kinetic 
approach A well defined statistics can be de- 

rived for any particular pair of values (M,N) indicating 
the degree of the polynomials, in the velocity variable, 
used to describe the drift J and diffusion D coefficients. 
The Tsallis classical and quantum statistics are related 
to the pair (0,1) (in this paper we do not derive the drift 
and diffusion coefficients of the solar core, this will be 
investigated elsewhere, rather we use straight way the 
generalized Tsallis statistics). 

In addition to this, let us recall that Clayton 
and collaborators J^,|22| suggested, on phenomenological 
grounds, that neutrino counting rate will be much re- 
duced if the high energy tail of the MB distribution of rel- 
ative energies is depleted, the depletion being described 
by the factor exp{— 5 (E/k B T) 2 } with a suggested value 
S » 0.01. 

The distribution that we have derived in Ref. |]l2| , |l3| , 
of the Tsallis type, can take into account this behavior, 
certainly due to the long-range gravitational interaction. 
The Tsallis parameter q can be related to the Clayton 
parameter <5 — (1 — (?)/2; our constraints impose a value 
of S slightly different from the value argued by Clayton. 
We obtain that the value of 5 must satisfy the range of 
variability 0.02 < S < 0.05. 

To show the analogies of our approach with the distri- 
bution suggested by Clayton, we introduce, by means of 



an expansion, an approximated expression of the Tsallis 
distribution which is correct in the range of energies of 
interest here. 

The Tsallis distribution function reduces to zero at E w 
(10 -7- 25) k B T depending on the value of 6; therefore, 
high energy collisions are greatly reduced or absent and 
the neutrino flux from 7 Be and 8 B reactions is smaller 
than foreseen by standard theories, the importance of the 
neutrino flux from pp reactions is consequently increased. 

In this work we show that the suggestion of Clayton 
and collaborators is well founded on theoretical grounds 
(long-range gravitational interaction), it is equivalent to 
using the distribution we have derived elsewhere (l^ [lj|] 
and that this distribution is well motivated within the 
Tsallis statistics. The neutrino counting rate, compared 
to the standard predictions, can be explained using the 
above prescriptions. A complete treatment of this sub- 
ject can be carried on only within a complete solar model 
code [p3|-p5|, it will be done and reported elsewhere. 



Following Tsallis |l|,|| Plastino and Plastino || and 
Boghosian jlTj , the distribution function / of stellar poly- 
tropes is 



/oc[l + (g-l)(a + /3e)] 1/ ^ , 



(3) 



where q is the Tsallis parameter and e = E-\-<$>(r), where 
E is the cm. kinetic energy, <E>(r) is the gravitational 
potential, a = —/3fi, where fi is the chemical potential, 
(3 = l/(k B T). The distribution (3) becomes the MB 
distribution exp[— (a + (ie)] as q goes to one. Let us 
define the relative gravitational potential W 



<3> = + $ 



(4) 



where $0 is a constant to be chosen such that \& vanishes 
at the edge of the system Q. The following relation 
between the relative potential \& and the density p holds 



,7-1 = Xzl* 



/C7 



(5) 



where JC = V/p 1 , V is the pressure and 7 is the adiabatic 
parameter. 

The polytropic index n is defined by 



7 = 1 + 



1 



(G) 



We introduce the relative energy £ = "J" — E, the distri- 
bution / can be written as 



/ oc [1 + (q - 1)(q + 0*0) - (<7 - l)[3£] ini ~ q) 



(7) 



The quantity a can be chosen, without loosing generality, 
in such a way that p7fl 



1 + (<Z-I)(a + /?$ ) = 



(8) 
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We want to compare the Tsallis statistics / to the dis- 
tribution introduced by Clayton et al. p2| which is equal 
to the MB distribution C exp(—(3E) times a correction 
factor 



f _ Ce~ f3E e v{E) = e Po-PiE-f3 2 E 2 



(9) 



where we set /?i = [3 + B (B is a constant). 

We observe that the Tsallis distribution can be written 

in a classical form / oc e~ c where 



1 



l-<? 



{log[l + (g-l)(a + /3e)]} 



(10) 



The function e may be interpreted as a generalized energy 
that takes into account many body collective interactions 

H- 

The logarithmic term in Eq.(10) can be expanded in pow- 
ers of (q—l)(a+(3e) with the conditions | (q— l)(a+/3e) |< 
1 and (q- l)(a + /3e) + -1. 
We impose that 



a + P$ = 



(11) 



to allow that /3i = (3, as required for physical reasons in 
the approximated distribution by Clayton ||^2| . 
Infact, B does not have any effect and can be ignored 
because we maintain the solar luminosity at its known 
value and the increase of the central temperature T to 
counteract the effect of B on the power of the sun raises 
the 8 B neutrino flux back to the value it had at B = 0. 
The special condition (11) or a = /3(*f? — <3?oX which is 
a constraint characteristic of the solar core, shows that 
a is a function of \& or (£ + E), therefore the compar- 
ison between two different but equivalent expressions of 
/ done to derive a relation between n and q in Ref.s 
P,p7| is not allowed in this case. Finally, we find that 
the parameters /3o, Pi, fii are: f3o = logC, (3\ = (3 and 
(3 2 = {l-q)^/2 = 5(3\ 
The above constraints (8) and (11) imply 



1 + (<? - l)/3# = , 
and from Eq.s (5), (6) and (12) we obtain 



T T 

or 5 (n + 1) 



n+ 1 



2 ' 



(12) 



(13) 



with t = k B T p/V. This relation, which links n and q, 
differs from the relation given by Plastino and Plastino 
H or from the one given by Boghosian (l7]] because of 
the constraint (11) valid for the solar core. 

In Ref. |(| it is shown that, if we select n — 3, the fol- 
lowing values of p, V and T of the solar core can be de- 
rived: the proton density p = 53.47 gr/cm 3 = 0.32 10~ 13 
protons/fm 3 , the pressure V = 0.77 10 -16 MeV/fm -3 
and the temperature k B T = 1.034 10 -3 MeV which is 
slightly lower than the temperature fixed by other models 
(1.29 10 -3 MeV). These last two figures are determined 



supposing a particular central composition; changing the 
composition by increasing the He concentration a greater 
temperature can be reached. 

With this selection of values we obtain r = 0.43 and 
8 = 0.05. 

In the solar core the value of the different physical 
quantities of interest in this work have been reported 
by Ichimaru S to be: k B T = 1.29 10~ 3 MeV, (3 = 
0.77 10 3 MeV^ p = 56.2 gr/cm 3 , which is the 36% of 
156 gr/cm 3 , for the free protons and V = 2.12 10~ 16 
MeV/fm 3 . 

By inserting these values into Eq.(13) we obtain (with 
n = 3) r = 0.2 and S(n +1) = 0.1; then we can fix 
the central composition to maintain the magnitude of r 
to be 0.2 and select different couples of values of n and 
<5, e.g.: n = 3/2 and 8 = 0.04, n = 3 and 8 = 0.025, 
ri = 5 and 8 = 0.017. In conclusion, we may expect that 
the appropriate value of 8 be in the range between 0.02 
and 0.05. The value proposed by Clayton (8 = 0.01) 
gives (when r = 0.2) n = 9 which is unphysical, because 
greater than 5. 

The approximated distribution function with 8 — 0.02 

is 



approx 



= c 



,-0.774 B ko 



(14) 



where C is the normalization constant. It is easy to ver- 
ify that the correct Tsallis distribution reduces to zero 
at E = k B T/(28) forbidding the existence of ions with 
energies greater than 10 4- 25 times the sun temperature. 
We wish to recall that the Eq.(14) represents a Druyven- 
stein distribution p7| . 

All quantities relevant in solar physics and solar neutrino 
emission can be expressed and evaluated as functions of 
the parameter 8. We report here, for instance, the ther- 
monuclear reaction rate corrected respect to the MB rate 
r MB by the depletion factor and calculated in Ref. (p2[ 



(15) 



where E is the most effective energy (E = 4.5 k B T for 
pp reactions) and A is a function of E , k B T and 8. 

The neutrino counting can be reduced sensibly above 
few k B T with some changes in the solar model param- 
eters. These can be compensated by small changes in 
the initial He concentration. Most of this reduction has 
come at the expenses of the 7 Be and 8 B neutrino fluxes. 
Counting rates within the solar model used by Clayton 
and collaborators can be extrapolated from the curves re- 
ported in Fig. 2 of their work J2^|. A value of 8 different 
from zero makes a star more luminous and reduces the 
rate of energy production at a given temperature. The 
solar core contracts to higher temperature. As shown by 
Clayton et al. the increase of temperature at given so- 
lar luminosity does not increase the neutrino fluxes that 
decrease with 8. 
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We do not discuss further, in this work, the measured 
results and the predictions of the neutrino fluxes; we 
leave complete and more definitive discussion after the 
proposed distribution will be tested within the available 
solar models |^,^,^5|,^8) . Of course, the results reported 
in this work depend on the values of the parameters of 
the solar core one takes as input. We can expect that 
the trend of definitive results will be on the line of the 
present description. We hope that the content of this 
work could be useful to the operating and proposed solar 
neutrino and underground nuclear astrophysics experi- 
ments HQ. 
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